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The spin evolution of a Bose-Einstein condensate starting from a mixture of two or three groups 
of 52 Cr (spin-3) atoms in an optical trap has been studied theoretically. The initial state is so 
chosen that the system does not distinguish up and down. In this choice, the deviation caused by 
the single-mode approximation is reduced. Moreover, since the particle number is given very small 
(N — 20), the deviation caused by the neglect of the long-range dipole force is also reduced. Making 
use of these two simplifications, a theoretical calculation beyond the mean field theory is performed. 
The numerical results are help to evaluate the unknown strength go. 

PACS numbers: 03.75.Mn, 03.75.Kk 



I. INTRODUCTION 

Bose-Einstein condensate (BEC) of atoms with 
nonzero spin has greatly attracted the interest of both 
the experimentalists and theorists [l], H, 0, 13, IE B 13 m 
recent years. Four years ago, the Bose gas of 52 Cr atoms 
with electronic spin J = 3 and nuclear spin 1 = was 
condensed successfully |8|. Experimentally, the atomic 
spins of 52 Cr were frozen in magnetic trap, but freed in 
optical trap For these spin-3 bosons, the interaction 
between two atoms is specified by the strengths gs, where 
S = 0, 2, 4 and 6 are the total spins of the pair. All gs 
except the one for S = have been determined experi- 
mentally To fully understand the characteristic of 
52 Cr and consider the use of this material in application, 
it is important to measure go with the help of theory. 
Among all the rich physics of BEC. one attractive phe- 
nomenon is the spin evolution [{J, EH EH EH EH Il5| . 
It was found that for 87 Rb and 23 Na, the evolution of 
the average populations of spin components sensitively 
depends on the strengths of interaction [l(| E3 • There- 
fore, the interaction can be determined (or confirmed) 
by observing the evolution. However, for 52 Cr, related 
experimental data and theoretical analysis are scarce. 

In this paper, we consider the spin evolution of a mix- 
ture of two or three groups of 52 Cr atoms in an optical 
trap. Each atom of a group has the same Z component 
of spin, /x. Magnetic fields are used in the preparation 
of these groups, but are cancelled during the whole evo- 
lution. What we are interesting in is the effect of the 
unknown strength go on the evolution. When the con- 
densate is very dilute and the temperature is very low, 
the single-mode approximation (SMA, namely, the spa- 
tial wave functions of the atoms distinct in fi are ap- 
proximately considered as the same) has been used by a 
number of authors to simplify the calculation. This ap- 
proximation has also been adopted here. However, even 
the above conditions are satisfied, the SMA might not 
be good as shown in [ll|. The validity of the SMA de- 
pends on the total magnetization of the system. In the 
following discussion, a special initial condition with total 
magnetization zero is chosen. Accordingly, the system 



does not distinguish up and down. Therefore, the devi- 
ation caused by the SMA is expected to be considerably 
reduced. 

It is well known that the system of 52 Cr contains 
the long-range dipole interaction which is considerably 
stronger than those of 87 Rb and 23 Na. The dipole inter- 
action is in general very weak. However, due to being 
long-range, the combined effect would be large if numer- 
ous atoms are involved. Besides, the effect would become 
larger and larger as the evolution goes on. Accordingly, 
we consider only the condensate with much fewer parti- 
cles and the early stage of evolution. For this very small 
condensate, the effect of dipole force is much smaller and 
therefore can be neglected. Furthermore, since the parti- 
cle number N under consideration is so small, the mean 
field theory might not work very well. Thus a method 
beyond the mean field theory is used in the follows. 



II. HAMILTONIAN AND THE EIGENSTATES 

For a 52 Cr condensate with N atoms, when the dipole 
force is neglected, the interaction between a pair of spin- 
3 bosons i and j is denoted as f/,j = OijS(ri —Tj), where 
Oij = goPo + g%Pi + 9iPi + 9ePe- And P to P 6 are 
the projection operators of the ^-channels. Based on the 
SMA, each boson has the same spatial state When 
integration over the spatial degrees of freedom has been 
performed, the Hamiltonian reads 

h = u » = / E °« ' / = / \^)\ 4 ^- (!) 



l<3 



1<J 



After the integration, only the spin degrees of freedom are 
left in the Hamiltonian. To diagonalize the Hamiltonian, 
we use the following Fock-states as the basis functions 



\a) = IWiWM), 



(2) 



where a indicates the set {-/Vg, • • • , N^}, is the num- 
ber of bosons with spin component n, and jl means — ji. 



There exist two restrictions on the as follows. 



N = A 3 + N2 + Ni + N + Ni + N 2 
M = -3A S - 2A2 - N- 1 +N 1 + 2N 2 - 



-A 3 
3A 3 



(3) 



2 



where M is the total magnetization (a constant). In 
the Fock space, by using the factional parentage coef- 



where C 3 '^ 3 ^ v ^ is the Clebsch-Gordan coefficient. In 

••• ,JV',-1,— ,N',—1,— 

the label o... ^-2 ■■■ " > the superscript denotes a 
revised set of {A^, ■ ■ ■ , A^} by reducing both N', and 
N' v , by 1, and the subscript denotes a revised set of 
{N 3l ■ ■ ■ , N 3 } by reducing A^ by 2. When the two revised 
sets are one-to-one identical, the label is 1, otherwise it 
is zero. 

When both N and M are given, the dimension of the 
matrix is finite. After the Hamiltonian is diagonalized, 
the j-th eigenenergy Ej and the corresponding eigenstate 
tjjj are obtained, ipj can be expanded by the basis func- 
tions (or vice versa) as follows. 

a j 

where the coefficient c J a = {o:\ipj} is real. These eigen- 
states are used in the following description of evolution. 



III. EVOLUTION OF POPULATION OF SPIN 
COMPONENTS 

The spin evolution begins when the three groups of 
atoms are mixed together. All the atoms of the first 
group have p, = 3, those of the second have p, = — 3 
and those of the third have [i = 0. Therefore, the ini- 
tial state is just a Fock-state \I) = |AT 3 , 0, 0, N , 0, 0, N 3 ), 
where N 3 , N 3 , and Nq are the number of atoms in the 
first to third groups, respectively. \I) can be expanded by 
the series of ijjj with the coefficients c'j. Thus the time- 
dependent solution of the Schrodinger equation ^(t) de- 



ficient [19( | , the matrix elements of Hamiltonian read 



(4) 



scribing the evolution reads 

V(t) = e- lHt ' h \l) = c^e- lE ^l n i, 3 . (6) 
3 

Since are known constants determined by \I), Ej and 
ipj are also known. So the evolution can be fully under- 
stood. 

From ^(t), we define the time-dependent population 
P^(t) which is the probability of an atom in spin compo- 
nent y, at t. It reads, 

P^(t) = -^*(t)l<a„l*(*)> = K + 0£(t) (7) 

where and a M are the creation operator and annihila- 
tion operator of an atom in /1, respectively. 

K = E(4) 2 E(^) 2 ^ (8) 

j a 

(t) = 2 £ cost ^-^ cM £ cU §- 0) 

3<j' a 

is the number of atoms in fi within \a). Equation ([7]) 
contains two terms. The first one B 1 is only determined 
by \I) and is time-independent and therefore, appears as 
a background of oscillation. The second one 0^(4) con- 
tains the time-dependent factor cos[(Ey — Ej)t/h] which 
implies an oscillation upon the above background. At 
the beginning of evolution (i.e., t = 0), ^(0) = \I) and 
P^(0) = NpN. Moreover, if the initial state has the 
symmetry = N^, the population would also be sym- 
metric where (t) = P^(t). This is because in this case, 
the Z axis can be reversed. 



\H\a) = {N^N{N^N[N^\fJ2O tJ \N- 3 N 2 N 1 N N 1 N 2 N 3 } 



i<j 



ii<v\n'<v' S=0,2,4,6 



-,N u -2, 



+ (1 - 5^1^)5^ 



JK>K>n»(n, - 1) , 



Jn^N'^N's-I) , 



+ (1 S^vjS^v - "—,N„-1,—,N V -1," 
/ ■■■.N'-1.---,N',-1, 

+(i - a„„)(i - s^Jn^n^n^j^:..^^ 
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In the following calculation, N = 20 is given. The 
initial states are given as N± 3 — (N — Nq)/2 and 
N± 2 — N± x = 0, where Nq is even and is ranged from 
to N. In this choice, \I) is uniquely determined by 
Nq. Obviously, the system has the up-down symmetry. 
And the total magnetization M is zero, which is a condi- 
tion in favor of the SMA. The 10 -8 meV, A and sec are 
used as units for energy, length and time, respectively. 
The strengths g2, <?4 and ge are taken from [2(|, namely, 
they are —3.88, 30.80 and 59.64, respectively, whereas g 
will be given at a number of testing values. The average 
density is given as / = 2 x 10~ n . This value is simply 
evaluated under a model that the density is uniform in- 
side a sphere. A slight deviation of / does not affect the 
following qualitative results. 
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A. The background 



FIG. 1: (Color online.) The background of spin evolution, 
B^, against go where (a) fx = 3, (b) fx = 2, (c) fx = 1 and 
(d) fi — 0, respectively. Five cases of initial states where 
iVo = 0, 2, 10, 18 and 20, respectively, are given and marked 
by the curves. 
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It is proved that the background of oscillation in spin- 
evolutions of sp in-1 condensates does not depend on the 
interaction [21| . However, the argument of that paper is 
based on the uniqueness of the Fock-state with a given 
N, M, and No- Obviously, the uniqueness holds no more 
for spin-3 systems. Therefore, the knowledge of interac- 
tion might generally be obtained by observing B 1 ^. This 
is shown in Fig. Q] where the dependency on the initial 
state and on g is revealed. It is also shown that B 1 ^ would 
depend on go rather weakly if go is positive. In this case, 
the structures of low-lying eigenstates would be domi- 
nated by S = 2 pairs (because only in this kind of pairs, 
the two atoms are mutually attracted). However, B 1 ^ 
would depend on go rather strongly if go is negative and 
close to <72- In this case, the structures of the eigenstates 
would vary sensitively with gp due to the competition of 
the S = 2 and S = pairs [22| . As a result, there is a 
domain of sensitivity. If the realistic go turns out to fall 
in this domain, it could be determined by observing B 1 ^. 



The background can be rewritten as B^ = J^j Wj Qj^ , 
where Wj = (c^) 2 is the weight of tpj in |7), = 
EM) 2N Z/ N = ^3K a Mj)/N is the probability of 
an atom in fx within rpj . Note that the curve with Nq = 
is much higher in Fig.[TJa), but much lower in Fig.[ljd). 
In this state, all the spins are either up (/i = 3) or down 
(p = —3) initially. Therefore, those ipj with a larger Q^3 
(i.e., having averagely more atoms for fi = ±3) would 
have a larger weight Wj . This fact explains why the 
curve with Nq = is the highest in Fig. [IJa), where 
the fi = 3 atoms are observed. Meanwhile, those ipj 
with a larger Qq would have a smaller weig ht Wj, which 
explains why the curve with Nq — is the lowest in 
Fig.Hd). 




0.002 0.004 



FIG. 2: (Color online.) The populations of spin evolution, 
P/,,(t), against t (in second) where (a) fx = 3, (b) fi — 2, 
(c) ft = 1 and (d) fx = 0, respectively. The initial state has 
-/Vp = 0. In every subfigure, the curves from the lowest to the 
highest have go from —40 to 40 with a step 10 as marked in 
(a). Each curve has been shifted upward by 0.2 more than its 
lower adjacent to guide the eyes. 



TABLE I: The first minimum and the second maximum of 
P/(t) denoted as -P^a an d f/.s, respectively, for No = 0. 
T/ a and T/ b are their locations (in ms). 
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1.84 


2.62 


4.88 
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1.94 
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B. The oscillation 0£(t) 

In Eq. (|5J), the time factor (Ej> — Ej)/h is in general 
not a multiple of integer among all pairs of j and j'. 
Therefore Plif) is non-periodic which is illustrated in 
Fig.H 

In the following discussion, we firstly focus on the case 
where TVg = 0. This case also implies that P± 3 (0) = 0.5 
at the beginning. Afterward, due to the appearance of 
other components, P± 3 (t) go down from the maximum 
0.5 as shown in Fig. [2^l, while others P^(t) go up from 0. 
And they fluctuate around the backgrounds. Let the first 
minimum of P^(t) in Fig. [DJa) be denoted as A located 
at t = A , and the second maximum in the figure be 
denoted as B located at t = B . Related data are given 
in Tab. HI It is clearly shown that from the locations of 
the maximum and minimum obtained via the theoretical 
calculation, the strength go can be determined once the 
realistic locations are experimentally measured. Addi- 
tional information can also be extracted from Figs. HJb) 
and[2jd). For example, the first peak of Ph(t) with fj, = 1 
or appearing in the earliest stage of evolution can help 
to discriminate g . 

For comparing, the evolutions for both Nq =10 and 20 
are illustrated in Fig. [3l It is shown that the evolutions 
are no more sensitive to go ■ Thus we conclude that Nq = 
as shown in Fig. [3 is a much better choice. 
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(c)N =20, n = 3 
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FIG. 3: (Color online.) The same as Fig. [2]expect that Nq = 
10 for (a) fi = 3 and (b) ju = 0, and Nq = 20 for (c) fi = 3 
and (d) \i — 0, respectively. 



IV. CONCLUSION 



We study the spin evolution starting from a mixture of 
two groups of 52 Cr atoms, which are fully polarized but in 
reverse directions and contains only a few particles. And 
we find an effective way for determining the strength ga- 
in this way, the deviations caused by the SMA and by 
the neglect of the dipole force are reduced. Accordingly, 
the theoretical approach becomes much simpler and a 
calculation beyond the mean field theory is performed. 
The numerical results show that the knowledge on g can 
be thereby extracted. Nonetheless, the above theoretical 
calculation can only provide a rough evaluation of go- For 
an accurate determination, more precise theory beyond 
the SMA and with the dipole force taking into account 
is necessary. This will lead to a great complexity, and 
hopefully can be realized in the near future. 
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